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ABSTRACT 


A nonstationary Markov model was developed relating weapon system 
effectiveness (WSE) to the time sequence of casualties observed in a 
two-sided, heterogeneous force, land combat simulation or field experi- 
ment. Parametric tests of hypothesis are used us analyse the relation- 
ships between WSE and engagement range and the numbers and types of 
combatants in the engagement, Maximum likelihood estimators of WSE 
demonstrated by a weapon system over the course of an engagement are ob- 
tained, Utilization of the Markovian model and WSE estimates in a low 
resolution simulation to investigate the impact of changes in force 


mix, speed of advance, and initial engagement range is discussed, 
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NOTATION 


A: denotes force A with c combatant types. 
A, : enobeswiweombatant type 1, 1=1,...,C. 


A(s,M,1;t,k)= the probability that the A combatants inflict k B casualties 
in the time interval feel, given M and N survivers at time s, 
A, (s,M,K;t,k)= the probability that the A, combatants inflict k B casual- 


ties in the time interval Cones given mM and N survivers at time s. 


A, (som, 15 t,k)= the probability that the A, combatants inflict k B, casu- 
alties in the time interval (AA given m and n survivers at time s. 
c d 
a(t,m,m)=a[ x(t),m,n ]= E Y a, ¿(t,m,5). 
i=1 3=1 79 
a, ¿Ctam,n)=a, [T(t), mn =the function used to define the probability of 
the As combatants inflicting Be casualties given m and n survivers at time 


t, denoted as an attrition function., 
B: denotes force B with d combatant types, 


B.: denotes B combatant type j, j=1,...,d. 


B(s,m,n;t,k)= the probability that the B combatants inflict k A casualties 
in the time interval (s,t|, given m and n survivers at time s. 


B (s,m,njt,k)= the probability that the B. combatants inflict k A casual- 


ties in the time interval ear given m and n survivers at time s, 


B (som mjt,k)= the probability that the B. combatants inflict k A. casu- 


alties in the time interval (Crue given Mm and TI survivers at time s. 
c d 
b(+,m,1)=b r(t) mn E Y b. (tm, i 
i=1 j=1 J 
b.,(t,mn)=b,.[r(t) msn] the function used to define the probability of 
the B. combatants inflicting As casualties given m and n survivers at time 
myecaenoted as an attrition function. 


K(1)= the number of casualties observed in the 1th engagement of a simu- 
lation or field experiment. 


L= the number of engagements or trails of a simulation or field experi- 
ment observed, 





L(p)= the likelihood function for an observed sequence of casualties and 
times between casualties with the set Pp of unknown parameters. 


M(t)= a cXi random vector LM, (1) ], the number of A survivers at tine t, 
with realization m. 

M,(t)= a random variable, the number of A, survivers at time t, i=1,...,C, 
with realization M, > 

m1” 
the kth casualty in the lth engagement, 


the cX1 vector of A survivers after the k-ist casualty and prior to 


N(t)= a dX1 random vector LN. (t) |, the number of B survivers at time t, 
with realization Nn. 


N.(t)= a random variable, the number of B, sürvivers ac tine Io- DER 


een realization a 


ng” the dX1 vector of B survivers after the k-1st casualty and prior to 


the kth casualty in the 1th engagement. 
P(s,m,nj;t,k)= the probability that the A and B combatants inflict a total - 
rcs on cach otber Tr the tine Tierna] (s,t ], given m and ü 
survivers at time s. 

p= the set of unknown attrition function parameters, 


T(t)= the matrix of ranges between opposing combatants at time t. 


as the matrix of ranges between opposing combatants at time t in the 


lth engagement. 

to,” the time the 1th engagement commences. 

tg” the time of occurrence of the kth casualty in the lth engagement. 
tol” he time the 1th engagement ends. 


1, if the kth casualty in the 1th engagement was a Er combatant 
SK me was inflicted by an As combatant, 

0; othewmise, 1=155 sac ae Jide 

1, if the kth casualty in the lth engagement was an A, combatant 
a jene was inflicted by a 2p combatant, 


0, otherwise, i=Í,...,;C and j=... sd. 





I. INTRODUCTION 


A. GENERAL 

The purpose of this paper was to develop a model and consequent 
procedures for the statistical analysis of land combat simulation and 
field experimentation data. The model was the result of an effort to 
provide a framework for the statistical analysis of data obtained from 
two-sided, real-time-casualty-assessment field experiments conducted by 
the U.S. Army Combat Developments Experimentation Command (CDEC). m 
ically these experiments chiefly addressed the measurement of weapon 
system effectiveness (WSE) and, specifically, the contribution of par- 
ticular weapon systems to total force effectiveness under a given sce- 
nario, In these experiments, WSE was measured in terms of the capability 
of a given weapon system or combatant type, rather than the capability of 
each individual combatant., 

The principal direction of this paper was the derivation and esti- 
mation of measures of WSE, recognizing that WSE is the result of inter- 
actions among a number -of factors describing the physical state of the 
engagement and the actions and decisions of the men employing the wea- 
pon systems, These factors were categorized as; the physical or per- 
formance characteristics of each combatant type, the numbers and types 
of combatants in the engagement, the engagement ranges, and other factors 
such as morale, terrain, and tactics, WSE is commonly assumed to be re- 
lated to weapon performance characteristics such as rate of fire and 
probability of a hit. Often denoted as measures of effectiveness, such 
factors may be experimentally determined, but their relation to WSE is 


not always clear. The capability of a given combatant type would also 
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be expected to depend on the number of weapons of that type engaged in 
the battle. Additionally, WSE may involve interactions with other wea- 
pon types and depend on the numbers and types of supporting and opposing 
weapons. Further, WSE may vary with engagement range through the depend- 
ance of certain weapon performance characteristics on range. Finally, 
WSE is related to various other factors, such as morale and tactics, 
which are difficult to quantify and are generally fixed by the simulation 
or field experiment scenario. 

Generally, loss exchange ratios (LER) have been the principal measure 
of WSE used in the analysis of the CDEC field experiments. LER are usu- 
ally defined as the ratio of enemy to friendly losses over the course of 
a battle or a given time period. LER are time average measures of rela- 
tive WSE and, hence, are sensitive to the duration of the engagement and 
the sequence of casualties experienced. Moreover, LER treat the rela- 
tions between WSE and the physical state of the engagement only implic- 
itly. Hence, LER may mask the possible dependance of WSE on such factors 
as force mixes, opening range of an engagement, or rate of advance. 

This study proposes quantifying WSE in terms of a En capability 
or potential for inflicting casualties. The proposed estimate of WSE 
explicitly evaluates the dependance of weapon system capability on range 
and the types and numbers of friendly and enemy combatants in the engage- 
ment as well as implicitly considering the relation between WSE and weapon 
performance characteristics. The estimate is an absolute measure of WSE 
as opposed to a relative measure as attained by LER, The WSE estimate 
complements rather than supplements the use of LER, as the former pro- 


vides an instantaneous measure of capability while the latter is a time 


average measure. 
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In addition to providing WSE estimates, the structure of the proposed 
model and the resulting estimates of WSE may be used to construct a lower- 
resolution, computer simulation of the engagement using minimal computer 
resources, The simulation may be used to investigate, within the general 
limitations of the initial scenario, the effects of varying such para- 
meters as force mix, initial engagement range, and speed of advance on 


the course of the battle. 























sensivity 
| analysis of | 
eee force mix, | 
etc, 


¡lower resol- $ 
on simul- {f 


[field exper- | 
iment or 
high 
resolution | 
p uain] 
B, THE MODEL 
The model is designed to investigate the attrition process in a two- 
sided engagement between heterogeneous forces, The attrition process is 
modeled as a nonstationary Markov process where the states of the process 
are defined as the numbers and types of surviving combatants and the 
ranges between opposing combatants. The transition probabilities between 
the Markov states depend upon the surviving combatants' capabilities or 
WSE which, in turn, depend on the performance characteristics of each wea- 
pon type, the numbers and types of surviving combatants, and the ranges 
between combatants. Thus, the model provides a means of relating the 
sequence of states and the times of state transitions observed during a 
battle to WSE, In turn, this allows for the estimation of WSE on the 
basis of the observed sequence of states and transition times. It is 
recognized that the course of the battle also depends on such factors as 


morale, terrain, and tactics: however, these factors are assumed constant 


over all replicates of the experiment and, as such, are implicitly included 
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in the transition probabilities. In this formulation, the estimate of 
WSE may be interpreted as the capability of a weapon system to inflict 
casualties and, specifically, as the rate at which a given weapon system 
inflicts casualties on a particular type of opposing weapon. 

In considering the validity of Markov assumptions, the nature of the 
combat process must be considered, An essential feature of any combat 
situation is its dichotomy. It consists of both an evolving physical 
system of weapons and their environment and a developing set of plans 
and decisions of opposing commanders. Given a state of the physical sys- 
tem and a particular set of plans and decisions, it may be argued that 
the evolution of the engagement depends only on the interactions of two 
opposing physical systems and is statistically deterministic--at least to 
the time of the next human decision. Of course the human element of plan- 
ning and decision making is anything but deterministic; under the same 
conditions, different commanders will often make different decisions. 
Thus a Markov model might prove a quite poor model of actual combat. How- 
ever, the use of decision rules in computer simulations and, to a lesser 
degree, the use of a scenario in field experimentation nullifies this 
human variability. In-the cases of simulation and field experimentation, 
the combat process has already been abstracted to a process with little 
or no significant human variability and whose future states depend only 
on the interactions of the opposing physical systems [ Koopman 1970 |. Ta 
addition, the Markov assumption implies that the future behavior of the 
engagement depends only on its current state, not on its past history. 
That is, if the course of the battle depends on the types and numbers of 
combatants, only the current numbers of survivers is important, not the 


sequence and times of casualties in reaching the current state. 
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The model is based on the standard assumptions for a multidimensional 
Markov death process: the independance of non-overlapping time intervals, 
the near zero probability of two or more casualties occurring simultane- 
ously, and, for all firer-target combinations, the existence of probabil- 
ities of a given firer type inflicting a casualty on a given target type. 
Given these assumptions, it is possible to develop, as a function of WSE, 
Chapman-Kolomorgov equations for the probability of observing the occur- 
rence of any particular vector of casualty types and numbers, Homer 
the resulting difference-differential equations are susceptible to a gen- 
eral solution only in special cases | Clark 1969 |. Fortunately, the gen- 
eral solution to the differential equation for the occurrence of zero 
casualties in a time interval is trivial. From this probability it is 
possible to obtain the general distribution of the time between casual- 
ties and formulate the likelihood function for an observed sequence of 
casualty types and intercasualty times. Consequently, maximum likeli- 
hood estimates of WSE or the rate at which a given weapon type inflicts 
casualties on a particular opposing weapon type may then be obtained from 


the likelihood function and simulation or field experimentation data. 


C. WSE ESTIMATION 

The quantification of WSE is a two-step procedure. Initially, the 
general nature E the functional relationship between WSE and the numbers 
and types of combatants and range must be determined. Secondly, given a 
general relation between WSE and the state of the system, the problem be- 
comes one of estimating the unknown parameters of this relationship. For 
example, if it is determined that the effectiveness of a tank platoon 


against a given antitank system is proportional to the number of surviving 
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tanks, n, and increases linearly as range, r, decreases, the general form 
of the tanks' WSE would be: WSE = n(a~br), where a and b are unknown 
parameters to be estimated, 

The general form of the functional relationship of WSE to the state 
of the system is obtained via parametric tests of hypothesis based on the 
asymptotic distribution of the generalized likelihood ratio. Essentially, 
in this step various general forms of WSE are postulated and the form 
which most closely correlates with the data is selected. Note that this 
procedure results in a statistical model of the attrition process, The 
general form of the estimated WSE equation is not necessarily the actual 
relation between WSE and the state of the battle; however, it is the form 
which provides the best available predictor of the course of the battle. 
The use of asymptotic or large sample distributional properties of the 
generalized likelihood ratio is justified even when relatively few ca- 
sualties of weapon type B, inflicted by weapon type A are observed since 
the analysis is also based on the observations of system states which led 
to casualties of another type or casualties inflicted by another weapon 


system type. 
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II. THE MODEL 


A. GENERAL 

A model was formulated to estimate weapon system effectiveness from 
the time sequence of casualties observed in a two-sided, heterogeneous 
force, land combat simulation or field experiment, It is implicitly as- 
sumed that each type of combatant Be quantifable capability to inflict 
casualties on each type of opposing combatant. This capability, which 
may be zero against some or all opposing combatant types, is denoted as 
an attrition function and is employed as a measure of WSE, The attri- 
tion function is assumed to be a function of weapon performance charac- 
teristics, the numbers and types of surviving combatants, and the ranges 
between combatants, The attrition function may reflect the values of 
other factors such as terrain, morale, and tactics, although these latter 
factors are assumed constant over all replicates of the engagement, 

The engagement attrition process was modeled as a nonstationary Markov 
process where the states of the process are defined by the numbers and 
types of surviving combatants and the ranges between combatants. The 
transition probabilities between the Markov states were assumed to de- 
pend on the surviving combatants' capabilities or WSE, which, in turn, 
depend on the performance characteristics of each weapon type, the num- 
bers and types of surviving combatants, and the ranges between combatants. 
With this formulation, a weapon system's attrition function or measure of 
. WSE represents the rate at which the weapon system is capable of inflict- 
ing casualties on a particular type of opposing weapon, 

The model is based on the standard assumptions for a multidimensional 


Markov death process: the independance of non-overlapping timeintervals, 
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the near zero probability of two or more casualties occurring simultaneously, 
and the existence, for each firer-target combination, of probabilities of 

a given firer type inflicting a casualty on a particular target type. Given 
these assumptions, Chapman-Kolomorgov equations may be developed, as func- 
tions of WSE, for the probabilities of observing the occurrence of any par- 
ticular vector of casualty types and numbers, The resulting difference- 
differential equations yield a general solution only in special cases, 
however, the general solution to the differential equation for the occur- 
rence of zero casualties in a time interval is straightforward. This 
solution is then used in the development of estimators and tests of hypo- 
thesis for the statistical analysis of simulation and field experimenta- 


tion data. 


B, THE ENGAGEMENT 

The two opposing forces are designated as force A and force B respec- 
tively. Force A has c combatant types with the ith combatant type des- 
ignated as Ass il... Co. Porce B has d combatant types with the jem 
combatant type designated as Da j=1,...,d. Combatant types are defined 
such that the separation between elements or individual combatants of a 
combatant type is negligible compared to the ranges to opposing combatant 
types. If the engagement scenario requires that identical types of com- 
batants eee widely separated groups, such groups are defined as 
distinct combatant types, A total of L independent replicates of the en- 
gagement are performed under the same scenario. 

Data available from each replicate of the engagement consists of the 
initial force levels, the type and time of occurrence of each casualty, 
and the type of combatant inflicting each casualty. In addition, unless 


any range dependence of the attrition process can be ignored, the relation 
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of the ranges between combatants to time must be, at least approximately, 
known. This relationship need not be constant over the replicates of the 
experiment and may be based on a known functional relationship, the sce- 


nario, or experimental data. 


CS ASSUMPTIONS 
Define the random variables: 


M, (t)= The number of A, survivors at time t, with realization n, , i=1, 
fee ,c, and 


N. BAY) The number of B. survivors at time t, with realization n;a ie 
Br ‚a. 


Define the random vectors: 


M(t)= The cX1 vector Dña (t) | of A survivors at time t with realization 
m, and 


N(t)= The dX1 vector In, (t)] of B survivors at time t with realization 
n. 


Define the variables: 
r, „(t)= The range between combatants A, and B, at tine t; 


T(t)= The matrix of ranges A ‚)] between opposing combatant types at 
Hime t, 


ae The maximum range at which combatant type A. can engage combatant 


type En and 
mai The maximum range at which combatant type B. can engage combatant 
type ae 

Define the probabilities: 


A, ,(s.,n5t,k)= The probability the A, combatants inflict k B. casual- 


ties in the time interval (s,t], given m and Ñ survivors at time s, and 


B ¿¿ (S/M, M3, k)= The probability that the B, combatants inflict k A, ca- 


sualties in the time interval (s,t], given f and m survivers at time s. 


Then the assumptions discussed above may be formally stated as: 


1. Events in non-overlapping intervals of time are independent. 
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2. For a given number of surviving combatants, casualties occur 
independently of each other in the infinitesimal time interval (t,t+h]. 
3. For a sufficiently small length of time h, there exist attrition 
functions a, ;LAt), MN] and bt), 1m,n] y i=l4... qe and j=l, e dice 
taat the probabilities of A. combatants inflicting k P casualties in the 
time interval (t,t+h] are given by: 
2. A, (tom nytth,t)=a, [r(t) msn Jh+0(n), 
ES 
D a EKOE 
k=2 7 
similarly the probabilities of ES combatants inflicting k A. casualties 
in the time interval (t,tth) are given by: 
C. B. (tmn; tth,1)=b. LE(t),m,n]h+0(h), 
Ms | 
Bee B..(t,n,n:eih,k)-0(h), 
k=2 J 
with the boundary conditions: 
A, ¿(t¡1,5;t+h,0)=1, if m,”0, n,-0, or 2, ¿(t) >r =: 


alJ 


B. (tmn; tth,0)=1, iL m=O, no or Ej) > 
A, (tm njtth,k)=0, if k>n.; and 


B (tm, nytth,k)=0, if k>m,. 


D. TRANSITION PROBABILITIES 
Let a, ‚(tım,n)=a, ‚Lr(t),mn] and b (tm,n)=b [FO mal. Define 
the probability: 


A, (s,m,n3t,k)= The probability that the A, combatants inflict k B ca- 


sualties in the time interval real given m and n survivors at time s. 


Then from the assumptions: 


d d 
A, (t,M,n;t+h,0)= II A, ‚(t,m,n;t+h,0)= II [1-a, .(t,n,n)ht0(h) ] 
p EM jean 


d 
zes a. „(t,m,n)_|hto(h), 
u 
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d d 
A. „(tm m,yn;tth,1)=> EN ¿ts ,n;t+h,1) II A, „(t ‚m,n;t+h,0) 
j= =1 * k=1 
kAj 
d 
> La, ¿Ct m,n)h+0(h)] 2 e: =a, „(tom ¿1 h+0(n) ] LE a, ¿ts m,n) | h+O(h), and 
j=1 
Wr 
N d 
ZA, (tom ‚n;tth,k)=0(h), N=Sn.. 
k=2 7 jai 4 


Analogous results follow for the attrition of A combatants by P combat- 


ants. Define the probability: 


A(s,m,n;t,k)= The probability that the A combatants inflict k B casual- 
ties in the time interval (s,t], given m and YN survivors at time s. 


Then from the jia above: 
A(tm,n;tth,0)= ir A, (t,m,1yt+h,0)= m ts a. „tm ‚n)_h+0(n)] 


i=1 i=1 jet > 


=1- Cs ; a, „(t,m,A) |nto(h), 
= 1 ja1 > 


A(+t,m,n;t+h,1)= Eh, (t,m,n;t+h,1) r A(t, m,n;t+h,0) 
i-i k=1 


kfi 
Cc d 
be 16 o a, (t,m,K) h+0(h) | 
í=1 j=1 me as 
ki 
LE y zas ¿(tm m mon), and 


i=1 -3=1 
N d 
Y A(t,m,n;tth,k)=0(h), N=E n.. 
k=2 j=1 9 
Again analogous results follow for the attrition of A combatants by B 


combatants, 


BE. CHAPMAN-KOLOMORGOV eo ERE 


c 
Define a(t,m,1)=x ; a, ¿ts m,ñ) and b(t,m,1) > > b.. (Emm anda 
al je =1 * i=1 j=l Ji 


the probability: 


no 





P(s,n,ü;t,k)= The probability A and B inflict a total of k casualties 
on each other in the time interval (Saale given M and n survivors at 
time sS., 

Then from the equations above: 
P(t,m,n;t+h,0)=4(t,1,15;t+h,0)B(t,m,n;t+h,0) 
=[1-a(t,m,n)h+0(9)] [1-d(+,1,1)h+0(h) | 


=1-[a(t,m,n)+b(t,m,1) Jpro(h), 


P(t,m,nytth,1)=A(t,m,n;tth,1)B(t,m,n;tth,O)+A(t,i,m;tth,O)B(t,m,m;tth,1) 
=[ a(t,m,n)h+0(h) ] [1-b(t,m,M)h+0(n) HLi1-alt,m,m)ht0(h) | [b(t,m,n)h+0(h) | 
-[ alt,m,n)+b(t,n,n) |h+0(h), and 
K E d 
E (t,m,n;yt+h,k)=0(h), K=E m,+ E n.. 
k=2 i=t * jet J 
Then the Chapman-Kolomorgov equation for the occurrence of zero casualties 
in the time interval (s,tth ], where h is an arbitrarily small length of 
time, is: 
P(s,m,n;t+h,0)=P(s,m,n;t,0)P(t,m,n;t+h,0) 
=P(s,m,n;t,0)[1-[a(t,m,m)+b(t,m,n) Jn+0(n) ] 
=P(s,m,n;t,0)-P(s,m,n;t,0)[[a(t,m,n)+b(t,m,1n) ]n+0(h) ], t>s. 
The Chapman-Kolomorgov equation yields the difference equation: 


P(s una syM,1Mjt,0) _ -P(s,1,n;t,0)[a(t,m,n)+b(t,m,n) |, 
and, when h approaches zero, the differential equation: 
<-P(s,m,13+,0)= -La(t,m,5)+0(t,m,1) P(s,1,1;t,0). 


Solving the above first order differential equation subject to the initial 


condition P(s,m,n;s,0)=1 results in: 


| t 
P(s,m,n;t,0)=exp_- f[a(y,m,n)+b(y,n,n) Jay]. 
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Similar differential equations may be obtained from the appropriate 
Chapman-Kolomorgov equations for the occurrence of any number of casual- 
ties in a time interval. However, the equations yield general results 
only with difficulty even for special forms of the attrition functions 
[ Clark 1969 ]. Fortunately, the probability of observing no casualties 
in a time interval yields sufficient information for the statistical 


analysis of simulation and field experimentation data, 


21 





III. STATISTICAL ANALYSIS 


A. GENERAL 

The estimation of WSE or weapon system attrition functions from 
simulation and field experimentation data was based on the maximum like- 
lihood criterion. The likelihood function for an observed sequence of 
casualties and intercasualty times was developed from the Markovian as- 
sumptions and Chapman-Kolomorgov equation by obtaining the distribution 
of the time between casualties and the conditional distribution of casu- 
alty type given the occurrence of a casualty. The estimation of WSE is 
then accomplished via a two-step procedure. Initially, the general nature 
of the functional relationship between a weapon system's attrition func- 
tion and the numbers and types of combatants and the ranges between com- 
batants must be determined, The form of the functional relationship of 
WSE to the state of the attrition process is selected by parametric tests 
of hypothesis based on the asymptotic distribution of the generalized 
likelihood ratio. Secondly, given the general form of the attrition 
function, the maximum likelihood criterion is used to estimate unknown 


parameters of the function. 


B. INTERCASUALTY TIME 
While the distribution of the number of survivors is relatively in- 
tractable, the distribution of the time between casualties can be obtained 
in a straightforward matter. Define the random variables: 
T,= The time of occurrence of the ith casualty, i > 1, with distribution 


function F(t), density £,(t), and realization t,; and 


C(s,t)= The number of casualties occurring in the time interval sten 


Ze 





Note that the event (T, >t) occurs if and only if the event [c(t, _,»t)=0] 
occurs, hence P(T, >t)=PLC(t, _,,t)=0]. If there are m and n survivors 
after the i-ist casualty, then 

1-F, (t)=P(T, > t)=PLC(t, ,»t)=0]=P(t, _,,M,M;t,0) 


1 
=exp| - J La(y,n,n)+b(y,m,15) Jay], and 
t 
1-1 


t 
St) = - PT, > t)=Lalt,m,n)+b(t,n,n) jexpL- J La(y,m,1n)+b(y,1,n) Jay]. 
il l 


C. CASUALTY TYPE 

If there are M and N survivors after the i-1st casualty and the ith 
casualty occurs at time as then, from the assumptions, the probability 
that the ith casualty was of combatant type ae and was inflicted by com- 


batant type As İs: 


a, .(t, .m,n)h o a, a(t, 0,7 
alt. ‚n,n)hrb(t. ‚n,n)h (t. ,m,n)+b(t. ,m,n) . 
Ñ í i i 


That is, given a casualty has occurred at time tio the random variable 
describing the type of the ith casualty and the type unit inflicting the 


ith casualty has a multinomial distribution. 


D. MAXIMUM LIKELIHOOD ESTIMATORS 
If there are ii and fi survivors after the i-ist casualty at time Te 
then the likelihood function associated with observing a ve Casualty 


inflicted by an A. combatant at time t is: 
ijti ns Br i as ae 
Gina ama, mn) t, m,n) Jex- [aly,m,n)+b(y,m,n) Jay] 
i t 
i-1 


ti 
p a lN La(y,m,n)+b(y,1,1) Jay]. 


i-1 
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To develop the likelihood function for an observed sequence of casualties 
and times between casualties, define the random variables and notation: 


1, if the kth casualty in the 1th engagement was a B. combatant 
Xj 47 and was inflicted by an A, combatant, J 


0, (Other see, 1=1,.5.,C 2nd Je, 


1, if the kth casualty in the 1th engagement was an A, combat- 
Y Sikl = ant and was inflicted by a Pe combatant, 


0, otheraze, I=ely.cn,ec and Jle oad; 
m7 The cXi vector of A survivors after ene k-1st and prior to the kth 
casualty in the 1th engagement, 
ny7 The dX1 vector of B survivors after the k-ist and prior to the kth 
casualty in the 1th engagement, 
= The time the 1th engagement commences, 
t,.= The time of occurrence of the kth casualty in the 1th engagement, 
t. y The time the 1th engagement terminates, | 
L= The number of independent replicates of the engagement observed, 


K(1)= The number of casualties observed in the 1th engagement, 


T, (t)= The matrix of ranges between combatant types at time t during 


engagement l, and 

p= The set of unknown attrition function parameters, 
At this point it is convenient to write the attrition functions as 
a, Lt, (+t) mi] and b [7 (4) mn] to reflect the possibility of changes 
in the relation between time and ranges over replicates of the engage- 
ment. If a sequence of casualties are observed over L independent rep- 
licates of the engagement, the associated likelihood function is given by: 

K(i) c d 


xX. e y ee 
= Il = — > JIRE 
1 = = a jee te) Mey pg J] ODL, Cty) sg hg IX 


ea EAE O) Ta Ta HE O) Ta T Jay he 
i) 
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o Lal ñ, (y) mol Hy Hole, (y) My Hay] E 
al 
The last term above is contributed by the time from the K(1)th casualty 
to the end of each engagement. For convenience and recognizing the loss 
of some information, the latter term is dropped in the remaining discus- 
sion. With this simplification, the natural logarithm of the likelihood 
ausetıon is: 
m E E EE galala, d n 
Bere) E- E Zx analna, IT (t-m an 4 
ee ee q a 
_ LK e A a 
Y 557 LD sy LE (ty ) M5, aad] 2 Be S Lal x, (y) m, 7, J F 
=1 k=1 t 
(k-1)1 
bZ (y) m, day] : 
Taking the derivative of In L(p)with respect to an unknown parameter 


gz which appears only in the attrition function aq. F(t) m1) and setting 


the derivative equal to zero results in: 


da [T(t),n,n] o I ka) x 


din L(p) din L(D pa k1 
dz = da rQ) mn] * dz 2 27 = Ce ieee ee 
pq sa l=1 k=1 “pa 1`”kl/’ kl!’ kl 


ACA E ua da LE, (y) Teee 
az l=1 k=1 t se 
(k-1)1 


dy=0, 


emanctLon of a [F(+) mm], Hence obtaining maximum likelihood esti- 
mators of parameters of a single attrition function involves working only 
with the attrition function of interest and is independent of the form 
and parameters of other attrition functions. 

Closed form solutions to the parameter estimators are attainable only 
in a few special cases and, in general, parameter estimation involves the 
solution of simultaneous nonlinear equations using numerical analysis 


techniques. However, the form of the likelihood equations readily lend 
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themselves to easy solution by simple analysis procedures. The parameter 
estimators for the various functional forms of the attrition function 
discussed in the appendices were found to rapidly converge to an appar- 
ently global solution using the basic numerical analysis technique of 
fixed point iteration {Conte 1972]. When an unknown parameter is common 
to several attrition functions, the derivative of In L(p) appears more 
formable, but does not greatly increase the difficulty of obtaining 
numerical solutions to parameter estimators. 

Confidence intervals or regions about attrition function parameters 
may be established via the large sample distribution of maximum likeli- 
hood estimates. The maximum likelihood estimates Py rDor Dy teem a 
sample size n of the parameters of an attrition function lm) are 
approximately distributed as a k- variate normal with mean vector 


(py +++.» Pp) and covariance matrix nR, where: 


nk -nEl In a(tsmin) 


mea 
AE 


a kXk matrix evaluated at p, =P, [Rao 1965]. 


E, TESTS OF HYPOTHESIS 

The general nature of the functional relationship between weapon sys- 
tems! Bi iiej on. functions and the numbers and types of combatants and 
ranges between combatants is SECS by parametric tests of hypothesis. 
This step of the analysis is based on hypothesizing various potential gen- 
eral forms of the attrition functions and then conducting sequential tests 
of hypothesis to select the form which most closely agrees with the ex- 
perimental data. The selection may be accomplished by simply hypothe- 
sizing a variety of functional forms and selecting from among them. Al- 


ternately, the selection process may start with a simple, say constant, 
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form of the attrition functions and sequentially test more and more 
complex forms until the addition of more parameters is no longer justi- 
fied by improvements in the fit of the data. In any case, this procedure 
results in a correlation model of the attrition process, 1.e., the form 
of the attrition function selected is not necessarily the actual rela- 
tion between WSE and the state of the battle; however, it is the form 
which provides the best available predictor of the course of the battle, 
The tests of hypothesis are based on the asymptotic distribution of 
the generalized likelihood ratio. Consider two hypothesized forms of 
the attrition functions with r and c unknown parameters respectively, 
mc, If L(P,,) is the likelihood function for the first set of attri- 
tion functions evaluated at the maximum likelihood estimates of the r 
unknown parameters and L(P,) is the likelihood function for the second 
set of attrition functions evaluated at the maximum likelihood estimates 
of the c unknown parameters, then -2 In L(p,)/1(p,.) | asymptotically ap- 
proaches a Chi-square distribution with r-c degrees of freedom [ Rao 1965]. 
In the comparison of two hypotheses where the numbers of unknown param- 
eters are equal, r=c, maximum likelihood criterion may be applied in the 
comparison of the two hypotheses. The hypothesis which results in the 
maximum value of the two likelihood functions evaluated at the maximum 
likelihood estimates of their parameters is selected. The use of asymp- 
totic or large sample distributional properties of the generalized like- 
lihood ratio is justified even when relatively few casualties of a given 
weapon type are observed since the test statistic is also based on the 
observations of system states which led to casualties of another type or 


casualties inflicted by another weapon system. 
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There are 2(cXd) attrition functions in the likelihood function and 
attempting to test composite hypotheses about this number of functions 
for the variety of permutations of possible functional forms would be 
time consuming at best. Fortunately, it is possible to reduce the test 
of hypothesis to a comparison of those attrition functions which have 
common parameters under one or both of the hypothesis, Consider the con- 
parison of two hypotheses which differ in the form of the attrition func- 
tion A pget(t) nn] and where the other functions are assumed to ie. the 
same, though perhaps unknown, forms with no parameters common to 
ang t(t) mn]. If the associated likelihood functions are L(p,) and 
L(p,.) » containing c and r maximum likelihood estimates of unknown para- 


meters of a lt) mn] respectively, then: 
-21n| LB) /L(B,,) J=21nL(P,)-21nL($,)= 


L A 


Zu 


Ii k-i X ngkılin a, a rn a a rel - 


r-c' 
“(KE DI l 


: 
l) Ea gA O) Ta arla IS 


and the test Statistic is ajtunetionsok a Ft) mn] only. If several 
attrition functions share common parameters, then Ye forms of these 
functions must be compared against any and the 
hypothesis test statistic can be reduced only so far as an expression in 


these functions. 
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IV, MARKOV SIMULATION 


Another consequence of the Markovian assumptions is the possibility 
of creating a computer simulation of the original simulation or field ex- 
periment using minimal computer resources, The Markov simulation is of 
lower resolution, but would be useful in investigating the effects of 
changes in those variables explicitly considered in the Markov model: 
the numbers and types of combatants, ranges, and the relationship between 
time and ranges. Potential applications includ: studies of the changes 
in engagement outcomes with changes in initial force mixes, initial en- 
gagement range, and speed of advance. As with all correlation cee 
caution must be exercised in considering parameter values outside the 
range of the data base; however, the simulation is intended more as a 
design tool for higher resolution experiments and should prove ante 
wen this role, 

The simulation is based on the Markovian assumptions of the model and 
the attrition functions estimated from simulation or field experimenta- 
tion data. An implication of the Markov model is that the battle may be 
regarded as a sequence of observations of two random variables: the time 
to the next casualty given a particular force level and the type of casu- 
alty observed given a casualty has occurred. The simulation simply in- 
volves specifying the initial forces, estimated attrition functions, re- 
lationship between time and ranges, and the engagement stopping conditions 
and then, using monte carlo procedures, generating a sequence of observa- 
tions of the Ko random variables discussed above. A flow chart of such 


a simulation is indicated in Figure 1, 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A, CONCLUSIONS 

The Markov model yields an estimate of a weapon system's effectiveness 
measured in terms of the rate at which the weapon system inflicts casual- 
ties on a given type of opposing weapon. The estimate is not a single 
measurement as are LER, but an evaluation of a continuous function which 
indicates the variation of WSE over the course of the battle as ranges 
and the types and numbers of combatants change. It is essential to dif- 
ferentiate between the estimated attrition functions and the values of 
these functions when evaluated over the observed values of r(t),m, and n. 
The analysis procedure guarantees only that it will select the variables 
and functional relationships which produce values of WSE which most 
closely agree with observed time sequence of casualties, not that it will 
select the attrition function which most closely agrees with the rela- 
tionship between WSE and the state of the attrition process. Essentially, 
the Markov model converts the information contained in an observed time 
sequence of casualties into a continuous function which, when evaluated 
over the observed sequence of casualties, provides a measure of the ca- 
pability demonstrated by a weapon system over the course of the engage- 
ment. For example, an analysis of 31 trials from CDEC experiment 43.6 
yielded the estimate of WSE for TOW equipped COBRA helicopters against 
aggressor tanks shown in Figure 2A, Effectiveness, measured as aggressor 
tanks destroyed per minute, is plotted against the range between the 


COBRA's and tanks and represents the average effectiveness observed over 


ehe 31 trials. 
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The estimates of WSE for individual weapon types may be combined to 
obtain the total effectiveness of a force against a given opposing weapon 
since, under the Markovian assumptions, the estimates are simply addita- 
tive. Given the validity of these assumptions, adding WSE does not in- 
flate the estimate of total force capability as the estimation procedure 
for individual weapon type WSE is based on this additative property. 
Figure 2B, again based on CDEC experiment 43.6 a plots the average 
observed total effectiveness of the defending force against the aggressor 
tanks. The plot also shows the contribution of the COBRA helicopters to 
the total force effectiveness. 

It should be emphasized again that the procedure discussed above does 
not predict WSE, but transforms information on WSE contained in an ob- 
served time sequence of casualties into a single dimension, a series of 
numbers which, in some sense, measure WSE and indicate how it varied over 
the course of a battle. An intermediate step in such an analysis is the 
development of a function which relates WSE to ranges and the numbers and 
types of combatants., Caution should be exercised in the use of this func- 
tion for predictive purposes as the Markov model guarantees only a corre- 
lation model of the observed engagement. For example, the function giving 
the effectiveness of the COBRA helicopters against aggressor tanks dis- 
cussed above indicates that WSE is proportional to the number of target 
tanks and increases linearly with decreasing range. As far as it goes, 
such a relationship is reasonable; however, the function suggests the 
COBRA force WSE is invarient as the number of COBRA's in the force changes. 
While possibly true, a more likely explanation is the fact that the data 
base used included very few COBRA losses and, to the Markov model, the 


number of COBRA's was essentially constant and, hence, not an important 
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COBRA Capability Against Aggressor Tanks vs Range 
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Variable in predicting the WSE of the COBRA force. Thus, using this 
function to predict the effectiveness of a different number of COBRA 
helicopters against tanks might be quite misleading; nevertheless, the 
estimates of WSE values for COBRA's obtained from the data are perfectly 
acceptable, This is simply an example of the risk, present in all cor- 
relation models, of attempting a prediction beyond the range of the data 
base. 

Another limitation of the Markov model, even when within or near the 
range of the data Dase, is the fact that the estimate of WSE effectiveness 
is a function of two random vectors, the vectors of surviving combatants, 
whose distribution is, in general, unknown. It is possible to make qual- 
itative predictions about the effects on engagement outcomes of changes 
in initial force mixes, speed of advance, and the like, However, quan- 
titative answers to such questions as expected losses or casualties in- 
flicted or the probability of winning can be addressed only by use of the 
Markov simulation. While the Markov simulation requires relatively little 
in the way of computer resources, its use is less desirable than an ana- 
INTE solution to these questions, 

Overall, the Markov model is primarily useful in quantifying WSE ob- 
served in simulation and field experimentation trials. If the limitations 
on its applicability are recognized, the model may be used to make predic- 
tions about the effects of scenario changes on engagement outcomes. Ap- 
pendices A and B contain examples of applications of the Markov model to 


the analysis of simulation and field experimentation data respectively. 


B. RECOMMENDATIONS 
Further work is needed in the validation of the assumptions used in 


the Markov model; however, given acceptance of these assumptions, the 
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model is ideally suited to the development of a standardized computer 
package for the quantification of WSE in the analysis of simulations and 
field experiments yielding time histories of casualty occurrence, Such 
an analysis package would use sequential tests of hypothesis to select 
from potential general forms of the attrition function and would develop 
maximum likelihood estimates of unknown parameters, The use of sequen- 
tial tests of hypothesis requires that the variety of potential general 
forms of the attrition functions be reduced to a computationally manage- 
able number, but remain comprehensive enough to include the variety of 
relationships between WSE and the state of the attrition process which 
experience indicates may exist. Appendix C includes an attempt to de- 
velop such a selection, The analysis capability should be complemented 
by the capability of using the estimatea attrition functions and casualty- 
time data to quantify and display the weapon capabilities indicated by the 
data. A crude version of such a package was constructed, using the at- 
trition function forms in Appendix C, for the analyses presented in Ap- 
pendices A and B, A similar standardized program could be developed for 
the Markov simulation for such applications as experimental design of 


higher resolution simulations and field experiments, 


5 





Appendix A. Analysis of Computer Simulation Data 


1. General 

This appendix presents the results of an analysis of computer simulation 
data using the Markov model, The data was generated using the Markov sin- 
ulation outlined above in order to exercise the analysis procedure in a 
situation where the attrition functions were known and the Markov assump- 


amos Valid. 


2o Simulation 

The simulation represented a hypothetical engagement between a 
stationary, defending A force and an attacking B force. The A force was 
composed of two combatant types: Ay and hos each with two elements. The 
B force was composed of two identical combatant types: B, and Bos each 
with four elements. The B, and B, combatant types were assumed identical 
to evaluate the consistency between the estimated measures of WSE for the 
two, The engagement opened at an initial range of 2200 meters with the 
B force closing at a constant velocity of 12 kilometers per hour to a 
range of 200 meters, where the engagement was terminated. The following 
hypothetical attrition functions were used for the simulation: 


a,,[7(t) mn =n Bil: i 
MEn] nn, 





a4 oL F(t) Suse] nl 40) , 


a, Lr(t) A n, ( ‚50-.00020r), 
A551 F(t) m,n |= nol .50-.00020r), 
b, ‚Lr(t) An n,m, ( 22000 EM 


d, pLE(t) 1,1 ]= (.30-.0001r), 





b,,Lr(t) TiS n,n, (.22-.0001r), and 
bool F(t) m,n |= (.30-.000ir), 
where r is the range between combatants in meters. 

Two independent sets of data were generated, the first with 20 repli- 
cates of the engagement and the second with 40 replicates to compare the 
Bis model under widely differing sample sizes. A total of 190 and 
390 casualties were observed in the first and second simulations respec- 
tively and resulted in the LER shown in Table 1, With 20 repetitions of 
the engagement, LER indicated that A, was more than twice as effective 


2 


against By as against B With 40 repetitions of the engagement, LER in- 


2 


dicated A, was nearly twice as effective against B, as against B In 


2 1° 
both cases, the hypothesized effectiveness against the two combatant 
types were identical and the anomalies in the LER were the result of the 


chance occurrence of disproportionate casualties during a part of the 


engagement, 


3. Analysis 

Figures three through eight compare the average values of the actual 
and estimated attrition functions for the observed time sequences of ca- 
sualties in the two simulations, The estimated values of WSE agree close- 
ly Ber the actual values and generally oa similar behavior over 
the course of the battle, Figure nine compares the average estimated 
values of WSE for A, against B, and B and B, and shows 


2 2 1 2 
that, as expected and unlike the LER, the estimated capability of each 


and A, against B 


weapon type is nearly the same against B, as against Boe 
Table two compares the hypothetical and estimated attrition functions 
for the two simulations. The correlation nature of the Markov model is 


apparent in that, despite the lack of significant differences in the 
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Simulation Loss Exchange Ratios 


time loss exchange ratio 
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Comparison of Average Observed Values of Actual and Estimated Attrition 
Functions, 20 Trials 
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Comparison of Average Observed Values of Actual and Estimated Attrition 
Mnet ons, 20 Trials 
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Comparison of Average Observed Values of Actual and Estimated Attrition 
Bmetions, 20 Trials 
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Comparison of Average Observed Values of Actual and Estimated Attrition 
Functions, 40 Trials 
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Comparison of Average Observed Values of Actual and Estimated Attrition 
monction, 40 Trials 
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Comparison of Average Observed Values of Actual and Estimated Attrition 
Functions, 40 Trials 
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Comparison of Actual and Estimated Attrition Functions 


Petr tion 
function 


a,,Lr(t) in] om 
a,oLr(t) mn] m 
ay, (F(t) sii] 
aso T(t) mn] 
b, ¿LE(t), m5) 
b, „LF(t) i] 
b„,Lrtt),m,n] 


ool T(t) siti | 


ris the range in meters between opposing combatants 


actual attrition 
function 


e= 34,10) 
n (+0) 
n, (. 50-.00020r) 
n, ( . 50- .00020r) 

m, (.22-.0001r) 
(.30-.0001r) 
m, (.22-.0001r) 


(.30-.0001r) 


estimated attrition functions 


20 trials 


n 
1 
mn, er (.578-,00022r) 


He )(.332) 
n, („914-.00040r) 
n(.47-.00020r) 
m, (.060) 
n„(.132-.00005r) 


m en ty) (151) 


Ir )(.196) 


Taole II 


W trials 


A 

ny (rra) (+280) 
i, m, )(.398) 
n,(.821-.00035r) 
aN e399 O 
n,(.548-.00047r) 
mo (+274-,00011r) 
m, (.097-.00003r) 


(.132) 





agreement of actual and estimated values of WSE between the two 
experiments, only the analysis based on forty repetitions was able to 
closely reproduce the attrition function forms and parameters used in 
the simulations. 

Given the validity of the Markov assumptions, this analysis indicated 
that even for relatively small sample sizes- approximately half of the 
sample size of CDEC experiment 43,6- the model yields satisfactory es- 
timates of the values of WSE demonstrated over the course of an engage- 
ment, However, considerably larger sample sizes are needed to develop 
a causative rather than correlation model of the attrition process for 


predictive purposes, 
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Appendix B. Analysis of Field Experimentation Data 


1. General 

This appendix presents the results of an analysis of CDEC experiment 
43,6 using the Markov model, Information and data on the experiment were 
extracted from the unclassified Experimentation Plan and Volumes I and II 
of the Final Report, USACDEC Experiment 43.6, Attack Helicopter-Daylight 


Defense, ACN 18171. 


2. Experiment 43.6 

CDEC experiment 43.6, phase II, was a field experiment simulating a 
two-sided engagement between an attacking tank heavy force with support- 
ing antihelicopter weapons and a defending force supported by antitank- 
missile firing helicopters. A primary purpose of the experiment was to 
evaluate the contribution of the missile firing helicopters to the cap- 
ability of the defending force. The nominal initial composition of each 


force is shown below: 


Defender, A force Attacker, B force 
combatant description number combatant description number 
type elements type elements 
A, COBRA helicop- 2 B, Tank 8 
ter with TOW 
A, Observation 1 B, REDEYE, SAM 1 
helicopter on APC 
À TOW on APC 2 B APC 2 
5 3 
Ay Tank 2 Bj 23MM SPAA 1 
De 57MM SPAA Í 
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The engagement commenced at a range of approximately 3750 meters with the 
attaching force closing with the stationary defending force at a velocity 


of about ten kilometers per hour. 


Dame Data 

The unclassified data available from experiment 43.6 consisted of 
casualty time histories from 46 repetitions of the above engagement per- 
formed in four similar areas. Data available from each experiment con- 
sisted of the initial force composition, the time and range at which each 
casualty occurred, and the type of and type unit inflicting each casualty. 
Since this analysis was intended as an exercise of the Markov model rather 
than a study of experiment 443.6, certain trials of the experiment were 
arbitrarily omitted from the analysis. Specifically, 15 trials, trial 
mompers 033, 184%, 194, 203, 214, 233, 254, 263, O11, 132, 161, 072, 082, 
171, and 201, were omitted due to instrumentation failure during the 
trial and/or the need for classified information to analyse the trial. 

The majority of the omitted trials could be included in the Markov analysis, 
but the purpose of this analysis did not warrent the effort. The omission 
of this data may have significantly biased the estimates of WSE, in that 
these trials nenne, included the ones most favorable to the attack- 
ing force, Again due to the limited one a this analysis, the data 
from different areas was pooled without regard to the effect of changes in 
terrain. 

One type of instrumentation failure which occurred occasionally in 
experiment 43.6 was the failure to render an assessed casualty inopera- 
tive which allowed this combatant to continue participating in the engage- 
ment. The Markov model can deal with such a failure directly due to the 


assumption of the independence of non-overlapping time intervals. In the 


49 





analysis of a trial containing such a failure, the casualty 1s assessed 
as usual, but the vectors of surviving combatants used in the analysis 
of the next casualty are left unchanged. 

Data casualty times were in clock, rather than lapsed time and, hence, 
the relationship between time and range was not directly comparable be- 
tween trials. The relationship between time and range was approximated 
by assuming the attacking force closed at a constant velocity of 10 kilo- 
meters per hour and least square criteria was used to convert casualty 
times into lapsed times from the beginning of each trial, where r(t=0)= 
3750 meters. Table 3 lists the adjusted time for the occurrence of the 


first casualty in each trial. 


4, Analysis 

The CDEC references on the experiment indicated that the combatants, 
except for Ay and Ba, were ineffective beyond a range of 2500 meters. 
However, the data records kills by other combatant types at ranges in 
excess of 2800 meters. The analysis is based on the assumption that all 
are ineffective at ranges beyond 2900 meters. 


combatants, except A, and B 


1 a 
Figures 10, 11, and 12A plot the average values of estimated WSE for each 
combatant against range aa time for the observed time sequences of casu- 
alties in the 31 trials. The data base did not include any casualties 
inflicted by a Ån, Bas Bus and Bg nor any casualties inflicted 
by A, on B, or B} on Ano Figure 12B indicates the contribution of com- 
batant Ay (COBRA) to the effectiveness of the A force against combatant 
B} (TANK). Tables four and five list the estimated attrition functions 


toe Lhe combatants, 
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Adjusted Time of Occurrence of First Casualty 


Times given in minutes 














experiment ty experiment ty 
023 7482 382 41,846 
123 9,696 102 An 
11 9.602 092 4,104 
154 5,646 all tt? 
173 5,434 021 8.538 
223 0.944 031 10.256 
24h 5,996 051 10.632 
083 3.428 062 0.924 
273 9.294 121 10.608 
293 11.394 1h2 3.620 
304 4.100 282 5.672 
314 12.070 212 4.322 
323 3.450 221 10.618 
090) 3.218 242 3.288 
34 - 5,354 261 10.154 
363 | 3.304 
Table III 


9l 





Average Observed Effectiveness of Combatant A, (COBRA) 


a, ¿LE(t) nn] 
US 


kil1s/minute ot] E 





<a A A PGI DARI EE IL ER DELLS PIE CL ETL TEN ER: : 
O 10 20 


time in minutes 
OO ee e e E a TO a aa 


range in meters 


Average Observed Effectiveness of Combatant Az (TOW) 
te 






TEG) nn] 
Alec yD 


kills/minute .3/ 


al 


time in minutes 
O OLA ne 21250 > AO 


range in meters 


Figure 10 


2 





Average Observed Effectiveness of Combatant A), (TANK) 


2, ¿_F(t) m1] 
JH1 100055 


ki11s/minute 





time in minutes 


range in meters 


Average Observed Effectiveness of Combatant B, (TANK) 


b,.[E(t) 1,7] 





i=1,3,4 
kills/minute 
b 
is? A 
b catar 
A a O 





a 1 20 
time in minutes 


ALTAS GRA 
range in meters | 





Figure 11 


I 





Average Observed Effectiveness of Combatant B_ (REDEYE) 
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Contribution of Combatant A, (COBRA) to Effectiveness of A Force Against 
Combatant B, (TANK) 
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All A combatant types appeared to distribute their fire uniformly, 


but the B, combatants seemed to distribute their fire only over the Ays 


3) and Ay combatants without engaging the A, conbatants, The geometric 


term or dependence on the number of targets which appears in attrition 


A 


fumetions By 4s a4, a s, and D indicates a high dependence of WSE 


33° 44 21 


cuate” target acquisition process, Attrition function Doo might be ex- 
pected to show a similar dependence had the data base included more than 
one A, (SCOUT) combatant, This result is consistent with the low alti- 
tude, pop-up tactics employed by the helicopters and the use of covered 
and/or concealed routes of approach by combatants B, (TANKS) and Ba (APC), 
The fact that this dependence did not appear in the other attrition func- 
tions may have resulted from the relatively fixed positions of combatants 
Az (TOW) and Ay (TANK) and a tendency for B, (REDEYE), By (23 MM), and Be 
(57 MM) to use more exposed moat ers in order to engage the helicopters, 
However, note that the data is insufficient to support such conclusions 
about the attrition of Bos Bus and Bs Since there is only one element of 
each of these combatant types. Similarly, it is impossible to determine 
the relationship between B, (REDEYE) effectiveness and the number of B, 


elements with a data base of one Bo combatant, Finally, due to the very 
small nunber of A, (COBRA) losses observed, it is doubtful whether the 


model would detect any dependence of A, WSE on the number of COBRA's, that 


dl 
is a44LE(t),m,11] might be mn, A, (.05-.00013r) rather than n,A,(.10-.00026r). 
For the most part, the attrition functions in Tables four and five pro- 
ibaply represent correlation rather than causative models of the attrition 
process and their use for predictive studies might prove risky. 

The plots of average observed WSE all tend to display a decrease in 


effectiveness toward the end of the engagement attributable to the dimin- 


ishing number of surviving firers and/or to the inereasimssscazelLy or 


a 





targets. Some of the plots show WSE to increase during the earlier part 
of the battle which is probably due to increasing weapon accuracy at 
shorter ranges and/or the increasing case of acquiring the remaining 
targets. This latter factor probably holds throughout the engagement, 
but is often outweighed by the diminishing numbers of firing and target 
combatants late in the engagement, | 

This discussion reveals some aspects of the comprehensive nature of 
the estimate of observed effectiveness of a combatant. This measure of 
WSE reflects not only the physical performance characteristics of a com- 
batant type, but also decreases in WSE due to the system's vulnerability 
to opposing weapons and consequent losses as the battle proceeds, More- 
over, the model reflects actual rather than potential WSE and may show 
WSE to decrease as the number of opposing combatants decrease due, say, 
to increasing difficulty in target acquisition and, of course, a system's 


WSE would go to zero when the opposing force was annihilated. 
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Appendix C 


1. General 

This appendix presents an attempt to develop a reasonably inclusive 
selection of potential general forms of the attrition function to con- 
sider in conducting sequential tests of hypothesis. The selection was 
designed to be sufficiently comprehensive to include known or suspected 
relationships between WSE and the attrition process, but limited enough 
to facilitate creating a standard program for the analysis of experimental 
data using the Markov model and sequential tests of hypothesis to select 
among the potential general forms of the attrition process. In addition, 
the consequent maximum likelihood estimators or simultaneous nonlinear 


equations which yield the estimators are indicated, 


Zeempeecneral Forms of the Attrition Function 

Consideration was limited to attrition functions with the general 
form of a product of a force level dependent term, f(m,m); a term repre- 
senting a combatant type's allocation of effort over potential target 
types, g(M,N); and a range dependent term, h[ x(t) |: 


a, LT (t) myn |-f(m,m)-X g(m,m) X hl x(t) |. 


Potential forms of the meee or force level dependent term were selected 
by analogy to existing Lanchester combat models. Forms considered were: 
a, Et) mA]=1 X g(m,m) x blx(t)], 
a, ¿LE(t) MM], X g(m,1m) x Hr()), 
a, sLt(t) mn j-n, X g(m,n) X hl x(t) ], and 


a, Ft), 12]=0m,n, X g(m,n) X nL x(t)]. 


2, 





Again by analogy to Lanchester theory, these forms of the arthi tiron 
functions are referred to as constant, square law, geometric law, and 
linear law forms respectively. 

The allocation of effort term was assumed to be either constant or 
proportional to the relative numbers of the opposing combatants, The 
constant allocation term simply represents the allocation of a constant 
amount of firer effort to each target type. The second form may be inter- 
preted as the allocation of effort in proportion to the relative numbers 
of each type of opposing combatant or as the result of a target acquisi- 
tion process where the probability of acquiring a given target type is 
proportional to the number of that target type. The attrition function 
may then be written as: 

a, ‚rm, nt xX 1X nis (0) or 
a, ¿LE(+),m, nl=f(m,n) X Gu, n y, > Mie o 


where Miss j=1, «es, d are unknown parameters to be estimated which re- 
present the relative importance of each opposing combatant in the allo- 
cation of effort, The last or range dependent term was assumed to be 
a linear function of range: h{ F(t) |= cee where a,b,C,»». 
are unknown parameters to be estimated, 

lt vas Belt that the er of attrition function forms resulting 
from the possible combinations of the terms presented above would prove 
computationally tractable and sufficiently general to describe most re- 
lationships between WSE and the state of the attrition process, Ap- 
pendices A and B contain examples of analyses conducted using the above 
forms of the attrition functions with h_Y(t) ] constrained to be of quad- 
tatic or lower Order de and ie in the allocation term one or zero, 


Natt s seg Oe 
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3. Maximum Likelihood Estimators 
The special form of the attrition function a, ¿LE(t) m,n jm; a is 


noteworthy in that the maximum likelihood estimator of the parameter a 


is given by: 
oe K(1) J L K(1) 
a ey ey ete St R 
et ket IKE peq ER DI 


which is the number of = casualties inflicted by A, combatants per time 
Breit of A, effort in inflicting these casualties, 
In general, for g(in,n)= 1 and hl x(t) |= a: 


Ki L k(1) 
= Een "a ran PR Ma) ha Hct 


a A | 
For g(m,m) =1 and h[T(t)]=(a'-b'r) = (a+bt), £ and b are given by the 
solutions to the simultaneous equationss 


Me Kil) a Se K(1) A 7 
a = > en 2 D f My (t -t NP and 
1=1 ket atbt, 4 ee N kl (k-1)1 


A E K(2) b x 277 ta J : Br 
a oe 


2 2 
ei, pies ee re ) 
2 ’ 
1= 2 1 kl 1l=1 k=1 kl KILO Kl oo e 


For g(m,n) = 1 and hl F(t) ] = (at-b'r-c'r*) = (atbtect”); a, b, and @ are 


given by the solutions to the simultaneous equations: 


Deere) ax. eek) 
A Iki 2 Q2 
a= 5 E aE / E E f(m., (t-t J 
nee Al kl’ kl’ kl “(k-1)1 
[ei ght 
L en Di L K(l) 
b= 5 ata; 5 4¢(m UE), ana 
ett ee NN a a ra) 
1=1 z 1 atbt, tote, 1=1 kei 
A ji S cx. Bet L i w C 
fer u en ym, (4) -t Ñ 
1=1 = atbt tote mee) K ar a (eine 
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